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ABSTRACT 

Transient microstructure in the diffuse interstellar medium has been observed towards 
galactic and extragalactic sources for decades, usually in lines of atoms and ions, 
and, more recently, in molecular lines. Evidently, there is a molecular component to 
the transient microstructure. In this paper we explore the chemistry that may arise 
in such microstructure. We use a PDR code to model the conditions of relatively 
high density, low temperature, very low visual extinction and very short elapsed time 
that are appropriate for these objects. We find that there is a well defined region of 
parameter space where detectable abundances of molecular species might be found. 
The best matching models are those where the interstellar microstructure is young 
(< 100 yr), small (- 100 AU), and dense (> lO" cm'^). 

Key vifords: ISM: structure - ISM: molecules. 



1 INTRODUCTION 

Observational evidence for microstructure on a solar- 
system scale in the diffuse interstellar medium has been 
accumulating for the last three decades. Measurements 
include VLBI of H atom 21 cm absorp t ion toward ex- 
tragalactic sources iRieter et alJ Il976t iDiamond et alJ 
Il989l : iDavis et alJIlQQgT abso rption of atomic lines toward 



transversely moving pulsars jFrail et al]|l99 6|). absorption 
differences toward stellar binaries jMeverfc Blades! 119961: 
iLauroesch et al ]j200tf) and toward multiple system s and 
clusters iLaur oesc h & Meyer 1999; Lange r et al] 
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in formaldehyde absorption towards extragalactic sources, 
the first report o f mole cular tracers of microstructure. 
iMoore fc MarscheJ lll995l) infer structure on scales of order 
10 AU, with number densities possibly larger than 10® 
cm~^, from further formaldehyde observations toward 
extragalactic sources. They also find spatial variations in 
formaldehyde and hydroxyl absorptio n line profiles towards 
the extended radi o gala xy 3C 111. iLucas fc Liszd lll993fl 
and lLiszt fc Luca3 ll2000h have used similar observations to 
explore the range of chemistry in diffuse and translucent 
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cloud s. It is also well known (IFalgarone. Pueet fc Peraul3 
I1992II that structure at larger scales (~ 0.01 pc) exists in 
nearby clouds. 

ICrawford] (jioO^) has recently detected absorption in 
lines of CH at the specific velocity (different to that of 
the surrounding medium) of a previously reported variable 
component in the line of sight to k Vel. Using a simple 
chemical model, Crawford obtained a value for the number 
density of the CH-containing microstructure toward k Vel 
that was consistent with the density inferred from measure- 
ments of the Cal/Call ratio at the same velocity component: 
riH ~ 10^ cm~^. The measured CH column density would 
require the path-length to be more than one order of mag- 
nitude larger than the assumed transverse size, supporting 
the view that the structure may be extended in the line of 
sight, either as a filament or a sheet. 

The most detailed st udy of mole cular microstructure 
to date has been made bv lRollinde et al. ( 2003 ). These au- 
thors made spectroscopic observations of the runaway red- 
dened star HD 34078 (whose velocity transverse to the line 
of sight is ~ 100 km s~^) over a 3-year period, supplemented 
by data from earlier epochs, to probe the foreground cloud 
distributions of CH, CH"*", CN and DIB carriers on scales 
from 1 - 150 AU. Their results show that the CH column 
density increased by 20% in 10 years while the CH"*" column 
density and two DIB strengths were unchanged. The CN 
column density shows a modest rise in this period. 

Although some authors have argued that some of the 
radio evidence for microstructure is flawed by poor data han- 
dling (e.g. Stanimirovic et al. 2003), the evidence in favour 
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Figure 1. Model geometry of microstructure at the velocity of 
the variable component in the line of sight to k Vel (different 
from that of the surrounding medium) with a depth transverse to 
the line of sight of 10 AU and an elongated length along the line 
of sight of 100 AU. The interstellar radiation field striking the 
microstructure transversely to the line of sight drives the chem- 
istry in the microstructure, its intensity falling off with Ay, trans- 
The contribution to column densities from species within the mi- 
crostructure is calculated along the line of sight through the struc- 
ture, over a visual extinction Ay^iong- 
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of the existence of microstructure now seems overwhelming. 
The general conclusion is that microstructure has a scale of 
10 - 100 AU, and has typical densities much larger than the 
ambient interstellar density. The transverse visual extinc- 
tion is therefore very small, and the gas in the structure can 
only be poorly shielded from the interstellar radiation field. 

The origin of the microstructure has been the su bject of 
intense study llPetv fc FakarondbOOOl: lGwinnll200ll) . Some 
authors argue that it cannot originate through non-magnetic 
hydrodynamics and it has been suggested that it may be 
excited in regions of high magn etic pressure by slow-mod e 
magnetosonic waves (Hartguis t. Falle &i WiLliamj l2003l) . 
Howe ver it is formed , it is overpressured and must be tran- 
sient. I^uin^^^^i] ll2003l) inferred from their observations 
that the variations were not the result of dense clumps 
passing through the line of sight. We are inclined to agree 
with their conclusion on dynamical grounds, and empha- 
sise that the transience of the microstructure implied by 
the high overpressure must affect the chemistry. It is un- 
clear whether these structures are extended, i.e. filamen- 
tary or sheet-like, or whether they are compact objects. Ex- 
tended objects could have longer path-lengths along the line 
of sight, if the orientation is suitable, and therefore the con- 
straints on density could be relaxed as compared to compact 
objects. The geometry may be related to the origin of these 
objects. 

The purpose of this paper is to explore the chemistry 
that may arise in the conditions apparently appropriate for 
the transient microstructure: rapid transition from relatively 
low to relatively high density, low temperature, very low 
transverse visual extinction, and subject to the normal in- 
terstellar radiation and particle fields. In Section 2 we de- 
scribe the model in more detail, and present in Section 3 our 
predictions of detectable chemistry as functions of time and 
depth into the slab. In Section 4 we describe the region of 
parameter space in which detectable molecular column den- 
sities might be found, and discuss the relevance of possible 
detections to the geometry of the structure. 



2 THE MODEL 

We model the region crudely as a semi-infinite slab illumi- 
nated from one side; the density rises as a step function at 



the edge of the slab. We compute the chemistry as a func- 
tion of time and transverse depth position within the slab, 
restricting our attention to the edge region (see Figure 0. 
What is unusual about this calculation compared to other 
studies of diffuse clouds is that the region of interest into the 
slab transversely to the line of sight is restricted to ~ 10 AU 
and to a very short evolutionary age, typically less than 10^ 
yr. These tiny spatial and time domains are usually ignored 
in studies of diffuse clouds (but note t he models proposed 
by ICecchi-Pestellini fc Dalgarnd teOOdl to account for H3" 
and other species in the line of sight towards Cygnus 0B2). 
We have therefore had to develop a chemical code in which 
particular attention is given to early time and edge effects. 
This code is a development of the pho todissociation region 
(PDR ) code written by Viti and Thi jPapadopoulos et alJ 

We investigate the chemistry in transient microstruc- 
ture using a PDR code that self-consistently determines the 
chemistry and thermal balance on a one-dimensional adap- 
tive spatial grid (Papadopoulos et al. 2002; to be described 
in detail in a subsequent paper). While the physical struc- 
ture of the slab is assumed fixed, the chemistry within it 
will be strongly time- and space-dependent. We have cho- 
sen gas phase elemental abundances, relative to wh, as fol- 
lows: He/H = 0.075, C/H = 1.79x10"'', N/H = 8.52xl0-^ 
O/H = 4.45x10"'', Na/H = 8.84x10"^, Mg/H = 
5.12x10"®, Si/H = 8.21x10"^, S/H = 1.43xl0"^ Cl/H = 
l.lxlO"^ Ca/H = 5.72xl0"^^ Fe/H = 6.19x10"^ The 
fraction of hydrogen initially contained in H2 is unknown; 
the chemistry will depend strongly on this factor. We note 
that the structures are dense and, if set up as proposed by 
Hartquist et al. (2003), would be found in post-shock re- 
gions. Hence, it is plausible to assume a significant fraction 
of H2 in the slab. We have arbitrarily set n(H2)/nH ~ 0.4, 
and we note that the chemical abundances predicted by 
the model should be proportional to this parameter. Since 
we are considering the additional contribution to molecu- 
lar species produced by chemistry in the microstructure, all 
other molecular abundances are set to zero initially. A total 
of 128 species are included in the model, connected through 
a network of over 1500 reactions, including ionisation and 
recombination of atoms. Freeze-out and mantle evaporation 
processes have been neglected due to the short timescales 
over which these structures evolve and their low visual ex- 
tinction. 

The filamentary or sheet-like nature of such objects is 
accounted for by assuming that the line of sight is within the 
slab and at a transverse depth of 10 AU from its surface (see 
Figure 0. Column densities are then simply proportional to 
the length of the line of sight within the slab. 

Since the physical conditions within these regions are 
uncertain, a large range of parameter space is examined to 
provide a thorough investigation of microstructure chem- 
istry. The effect of varying the slab density and age, as well 
as the environmental parameters of incident radiation field 
strength and cosmic ray ionisation rate ((") are all consid- 
ered, producing a four-dimensional grid of test parameters. 

Observational estimates of microstructure density gen- 
erally vary between 10'^ 5S riH ^ 10^ cm"'^, with some indi- 
cations of even higher densities. We have explored the model 
chemistry for this density range, and have made some ad- 
ditional calculations for the higher densities suggested by 
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Figure 2. Contours of log CO column density as a function of 
density and time for a microstructure of transverse depth 10 AU, 
elongated by a factor of 10 along the line of sight, and a cosmic 
ray ionisation rate of f = 1.3x10"^'^ s~^. From top to bottom, 
the three figures are for incident radiation field strengths of 1/3, 
1 and 3 Habing, respectively. 



Figure 3. Contours of log CH column density as a function of 
density and time for a microstructure of transverse depth 10 AU, 
elongated by a factor of 10 along the line of sight, and a cosmic 
ray ionisation rate of f = 1.3x10"^'^ s~^. From top to bottom, 
the three figures are for incident radiation field strengths of 1/3, 
1 and 3 Habing, respectively. 
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the formaldehyde observations of lMoore fc MarscherHl995l) . 
Slab ages of 1 - 1000 yr are considered. Radiation strengths 
and cosmic ray ionisation rates of 1/3, 1 and 3 times the 
standard interstellar values (1 Habing and 1.3x10^^^ s^^, 
respectively) are also considered. 



3 MODEL RESULTS 

The purpose of our study is to investigate whether there ex- 
ists a parameter space of densities, radiation fields, cosmic 
ray ionisation rates and dimensions where molecules can be 
formed in the variable components, and also to determine 
the time-dependent nature of such structures. In the next 
two subsections we make predictions, based on our models. 
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Table 1. Column densities (cm~^) for molecular species along the line of sight within a microstructure of transverse depth 10 AU 
and elongation factor 10 along the line of sight. The region is subject to standard interstellar radiation and particle fields of 1 Habing 
(1.6x10"'^ erg cm~^ and ^ = l.SxlO"^'^ s~^. Values in bold indicate potentially detectable abundances (with column densities 

> 10" cm-2). 
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Table 2. Column densities (cm ^) for neutral and singly ionised calcium along the line of sight within a microstructure of transverse 
depth 10 AU and elongation factor 10 along the line of sight. The region is subject to standard interstellar radiation and particle fields. 
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of the physical conditions and time-dependence of the mi- 
crostructure, respectively. 

In all the results reported here, we have assumed an 
elongation factor of 10; i.e. the path-length along the line 
of sight through the microstructure is 10 x the transverse 
depth, 10 AU. 

3.1 Physical Characteristics of the Microstructure 

The results of our chemical model suggest that only a few 
molecular species are formed in large enough quantities in 
the proposed scenario to be potentially detectable. They 
include CH, CHa, CO, Ca and OH. In our model, these 



species are only produced in great enough quantities to be 
detectable under certain physical conditions (see Tables 
and 121 and Figures |21 and |3J • In Table Q we list the column 
densities of CO, CH, Ca and OH as a function of density and 
time at a transverse depth of 10 AU into the microstructure 
slab, corresponding to a transverse visual extinction of ap- 
proximately 10~" to lO"'^ mags for number densities in the 
range 10^ to 10^ cm^"*, making the assumption that the 
elongation factor along the line of sight (described above) 
is 10. Column densities larger than lO'^^ cm~'^ are shown 
in bold in Table indicating a reasonable lower limit for 
molecular species that are potentially detectable. Assuming 
that the microstructure has a transverse spatial dimension 
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Figure 4. Fractional abundances of CO (solid line), CH (dotted), 
C2 (dashed) and OH (dash-dot) at a transverse visual extinction 
of Ay^trans = 10~^ mags into the microstructure, subject to an 
incident radiation field strength of 1 Habing and a cosmic ray 
ionisation rate of = 1.3xl0~^^ s~^. From top to bottom, the 
three figures are for microstructure densities of 10^, lO"* and 10''' 
cm~^, respectively. 
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of 10 AU, a minimum density of 10* cm~^ is necessary for 
most species to be detectable. Table |5| lists the column den- 
sities for Cal and Call. Since atoms have stronger oscilla- 
tor strengths than molecules, these species are detectable at 
lower column densities than the limit of 10^^ cm~^ imposed 
in Table Q Figures |21 and |3 show how the column densities 
of CO and CH behave as a function of density, time and 
radiation field. We find that CO is particularly sensitive to 



variations in the number density, and its column density can 
vary by over 2 orders of magnitude in the number density 
range considered here. Its sensitivity to the radiation field 
is more subtle: at early times (< 20 yr) and low densities 
(< 10* cm~'^), a weak radiation field (< 1 Habing) yields a 
longer timescale (~ 100 yr) for the production of detectable 
CO. CH, on the other hand, shows a strong sensitivity to 
the radiation field employed. In fact, if the field is stronger 
than 1 Habing, CH can only remain detectable at high den- 
sities (> 5x10* cm~^). In general, CH is formed and then 
destroyed very quickly. These trends restrict the parameter 
space for which both CO and CH can be present in the mi- 
crostructure at the same time. None of the species seem to 
have a strong dependence on the rate of cosmic ray ionisa- 
tion employed. 

The most remarkable result is that, even at the earliest 
times shown, if the number density is sufficiently high then 
potentially detectable column densities of the species listed 
above are attained. Thus, the high density drives a chemistry 
fast enough to overcome the losses due to the strong radia- 
tion field. Since the transverse visual extinctions considered 
are close to zero, the region is essentially unshielded. Note, 
e.g., that the to tal extinction alo ng the line of sight toward k 
Vel observed bv lCrawfordI (|20o3) is smaU {E{B-V) = 0.11) 
and that the transverse extinction at the variable component 
must be significantly less than this. By assuming that the 
transverse extinction is close to zero, we are therefore es- 
timating the lower limit of the additional contribution to 
column densities from species generated in the microstruc- 
ture. 

The adopted column density of H2 is, however, suffi- 
ciently large to ensure that the H2 photodissociation rate 
is too slow to cause significant loss of H2 on the timescales 
considered here. 

The fact that high densities are necessary in or- 
der to produce sufficiently large molecular abundances in 
the microstructure is co nsistent with the conclusions of 
iMoore fc MarscheJ (^9^ who find that, if the molecular 
variability is indeed due to density structure, the number 
densities need to be at least 10^ cm~^. I n fact, Marscher et 
al. (1993) and lMoore fc Marscheil |l993) detect both H2CO 
and OH in their time variability study. In our study, OH is 
present as long as the density is > 5x10^ cm~^ and does 
indeed show time variability (see Table 0. H2CO, how- 
ever, never reaches a high enough abundance to be de- 
tectable for the number densities considered here. In or- 
der to see whether densities as high as the ones consid- 
ered bv lMoore fc Marscher! ^1995^ could produce significant 
amounts of H2CO, we ran two more models at 10^ and 10^ 
cm""^ under standard interstellar conditions. We find that 
at 10'' cm~'^ H2CO is only just barely close to the lower 
limit of detectability for times ^ 500 yr, while at lO'' cm~'^, 
H2CO becomes det ectable after only 50 years. These models 
further support the lMoore fc Marscheil (^9^ findings. 

3.2 Time-Dependence of the Chemistry 

Tables and |5| and Figure 3] show how the different species 
behave as a function of time. We now briefly discuss the time 
variability of each species for number densities ^ 10* cm"''. 
CO takes ~ 5 - 10 years to form in significant amounts, and 
its abundance increases steadily with time. For CO to be 
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abundant, therefore, the microstructure would need to exist 
for at least 10 years. CH, on the other hand, forms in less 
than 1 year, but is subsequently destroyed (in part to form 
CO) within the first 100 years. C2 takes over 10 years to 
form in detectable quantities, then survives in the gas for 
at least 300 years, although its abundance starts to decline. 
OH behaves in a similar manner to CH in that it forms 
within 1 year but starts being converted into other species 
almost immediately and by 500 yr its abundance is too low 
to be detect a ble. C aH is present, confirming the detection 
of ICrawfordI i2002l) at the velocity of the microstructure. 
However, it cannot be used as a diagnostic test of the mi- 
crostructure age because its abundance is invariant. Cal, on 
the other hand, shows considerable time variability, though 
whether the variations are large enough to be detectable is 
less certain. Of course, these variations are not reflected in 
the column densities of CaH, since its abundance is several 
orders of magnitude higher. The CO column density un- 
der appropriate parameter choices generally exceeds those 
of other molecular species. Thus, for a number density of 
10* cm~^, by 100 yr the CO column density has built up 
to 2.6x10" cm-^ while those for CH and C2 are 5.8xl0"' 
and 4.0x10^"^ cm"'^ respectively. Although the amounts by 
which CH, OH, C2, etc. decline with time are very small 
and may therefore seem insignificant, they are indeed real 
and consistent with the gas phase chemistry that is taking 
place. The behaviour described above is expected for the 
usual gas phase chemical networks: CO and C2 column den- 
sities increase monotonically over the time span considered 
here. A feed molecule for the networks producing CO and C2 
is CH, and its column density declines slightly over this pe- 
riod, showing that CH is being converted into other species; 
no significant loss of CH by photodissociation occurs in the 
time interval explored here. 

In summary, for CO and C2 to form, the microstructure 
needs to survive at least 10 years. On the other hand, it 
cannot be much older than 50 yr if CH is indeed observed. 

Unfortunately, detection of CO by absorption in its res- 
onance line requires a UV spectrometer in an orbital ob- 
servatory; none is currently available. However, CH and C2 
are detectable from the ground. It appears from this work 
that these molecules should have comparable abundances in 
transient microstructure of number density ~ 10* cm^'^. 



4 COMPARISONS WITH OBSERVATIONS 
AND CONCLUSIONS 

This study was prompted i n par t by the molecular detec- 
tion obtained bv ICrawfordI (l2002h at the velocity (different 
to that of the surrounding medium) of previously reported 
microstructure in the line of sight toward k Vel. Hence we 
now briefl y compare our m odels with such observations. In 
Table 2 of ICrawfordI l)20n2h . the observed column densities 
of Cal, Call (two velocity components) and CH are listed 
as, respectively, 3.2x10", 2.0xl0" and 6.2x10^", 2.7xl0" 
cm~^. If we assume an elongation factor of 10 along the line 
of sight and a gas density of the order of ~ 10^ cm~^, the 
observed CH column density can be easily reached at very 
early times (~ 10 yr) and maintained. If instead the density 
of the microstructure is closer to 10* cm~^ then we infer a 
maximum age for the microstructure of about 50 yr, after 



which most CH would be destroyed. Hence observations of 
CH indicate that such filamentary or sheet-like structures 
are likely to be younger than 50 yr. The observed Cal and 
Call column densities can be produced in the model if the 
microstructure has a number density of ~ 5x10* cm~'^ and 
is older than 10 yr (see Table I^J. 

In order to further test the predicted molecular abun- 
dances, we conducted a high-resolution search for interstel- 
lar C2 at the velocity of the variable component toward k 
Vel using the Ultra-High-Resolution Facility (UHRF; Diego 
et al. 1995) at the AAT in February 2003. Sixteen indi- 
vidual exposures were made of the region containing the 
Q(4) line at 8763.751 A, with a total integration time of 4.9 
hours. Data reduction and calibration was as described by 
ICrawfordI il997l . 120021) . The final continuum signal-to-noise 
ratio was 480, but the Q(4) line was not detected. The re- 
sulting 3(7 upper limit to the equivalent width of the Q(4) 
line was 0.09 mA, yielding an upper limit to the J = 4 
column density of N{J = 4) < 1.95x10^^ cm~^ (using the 
oscillator strength adopted by Crawford 1997). At the high 
densities (jih = 10^ to 10^ cm~^) expected for these vari- 
able components, the rotational populations will be approx- 
imately thermal, and for T ~ 100 K we expect the total 
C2 column density to be ^(02) ~ 4 x N{J = 4) (e.g. van 
Dishoeck 1984). Thus, our observations imply ^(02) < 10^^ 
cm~'^ in the variable interstellar component toward k Vel. 
This is consistent with the model results given in Table 
for TVh ^10^ cm""^. More sensitive searches for C2 would 
appear to be worthwhile, as these would be able to constrain 
the model results at lower densities. 

The best matching models for observations at the ve- 
locity of the variable component in the k Vel line of sight 
seem therefore to be those for a young (~ 50 yr), transient 
(< 100 yr), high density (> 10* cm"^), small (~ 10 AU), 
slightly elongated (by about a factor of 10 along the line of 
sight) structure, immersed in an ambient radiation field. 

The main result of this paper is that essentially un- 
shielded interstellar gas can develop detectable column den- 
sities of molecules on very short time-scales (< 100 yr) for 
path- lengths on the order of 100 AU along the line of sight, 
if the density is ~ 10* cm"'' or higher. This time and den- 
sity domain is a region of parameter space not previously 
explored. The calculations presented here are pseudo time- 
dependent, in that the slab structure is fixed. 

If the cause of such microstructure is indeed magne- 
tohydrodynamic, then the filamentary objects likely to be 
responsible for the effects observed toward k Vel and else- 
where are probably waves moving at magnetosonic sound 
speed. More complex and r ealistic models are now being ex- 
plored dCarrod et al.ll2004^ . 
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